Intrinsic Spin Lifetime of Conduction Electrons in Germanium 
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We investigate the intrinsic spin relaxation of conduction electrons in germanium due to electron- 
phonon scattering. We derive intravalley and intervalley spin-flip matrix elements for a general spin 
orientation and quantify the resulting anisotropy in spin relaxation. The form of the intravalley spin- 
flip matrix element is derived from the eigenstates of a compact spin-dependent k-p Hamiltonian 
in the vicinity of the L point (where thermal electrons are populated in Ge). Spin lifetimes from 
analytical integrations of the intravalley and intervalley matrix elements show excellent agreement 
with independent results from elaborate numerical methods. 

PACS numbers: 85.75.-d, 78.60.Fi, 71.70.Ej 



OO 



o 

CZ5 



C3 



c 

o 
o 



> 

(N 
00 

in 
o 

(N 



X 



I. INTRODUCTION 

Group IV semiconductors are natural material choices 
for quantum and classical spintronic devices^— Hyper- 
fine interactions are suppressed due to the natural abun- 
dance of zero-spin nuclear isotopes. As a result, localized 
electrons have exceedingly long coherence times at low 
temperatures^— As for conduction electrons, space in- 
version symmetry precludes their spin relaxation by the 
Dyakonov-Perel mechanism^ The intrinsic spin lifetime 
is therefore relatively long, reaching ~10 ns at room- 
temperature in non-degenerate n-type silicon ,2r2£ Com- 
bined with the fact that silicon is the material of choice 
in the semiconductor industry, there is a wide interest in 
related spin injection experiments J^r— 

The motivation for studying spin injection in Ge is 
similar to Si due to their shared properties and the com- 
patibility with Si-based CMOS technology. Electrical 
spin injection and extraction in Ge have been recently 
investigated in lateral spin-transport devices with vari- 
ous doping profiles using nonlocal 19 and local22r— Hanlc 
measurements, as well as in heterostructure and nanos- 
tructure devices i 26 ' 27 Similar to direct band-gap semicon- 
ductors, optical orientation is an additional viable tool 
to investigate spin properties of electrons and holes in 
Gei^r— Unlike silicon, optical orientation in Ge is effi- 
cient because of the energy proximity between the direct 
and indirect gaps. Spin-polarized electrons are first pho- 
toexcited to the T valley and then they relax via ultrafast 
spin conserving scattering to the conduction band edges 
in one of the four L valleys (located ~140 meV below the 
zone center T- valley) <21 

Theoretical efforts in the early day o 35 ' 36 were moti- 
vated by low-temperature electron spin resonance ex- 
periments that studied the g factor and spin-lattice re- 
laxation of localized electrons in donor states) 3 - 7 - - — On 
the other hand, little attention was paid to conduc- 
tion electrons whose spin relaxation is mediated by the 
Elliott- Yafet mechanismi 40 ' 41 By analyzing the space in- 
version and time reversal symmetries of the L point, 
Yafet deduced a T 7 / 2 temperature dependence of the 
spin relaxation rate due to intravalley electron scattering 



with acoustic phonons^ 1 - Kalashnikov extended Yafet's 
theory to various statistical distributions and scatter- 
ing mechanisms.— Chazalviel investigated spin flips due 
to electron-impurity scattering using effective spin-orbit 
coupling parameters that resemble the treatment in III-V 
semiconductors M~ Most recently, Tang et al. have used a 
tight-binding model to calculate the intrinsic spin relax- 
ation of conduction electrons in strained Ge^ 

In this paper, we present a theory of spin-flip processes 
due to electron-phonon scattering in Ge. Two distinctive 
contributions are present in this work. First, we find 
the spin orientation dependence of spin-flip matrix ele- 
ments. This dependence leads to anisotropy in spin relax- 
ation and it is instrumental in analyzing measurements 
where the orientation of injected spins is set by the shape 
and magnetocrystalline anisotropy of ferromagnetic con- 
tacts or by the propagation and helicity of a circularly- 
polarized light beam. An interesting result of the analysis 
is that most of the intrinsic spin relaxation of conduction 
electrons in Ge can be explained by coupling of the lowest 
conduction band to the upper conduction bands (rather 
than to the upper valence bands which is the typical case 
in most semiconductors). The second contribution of this 
work is the derivation of a spin-dependent k-p Hamil- 
tonian in the vicinity of the L point (conduction band 
edge). This compact Hamiltonian model exquisitely cap- 
tures the signature of spin-orbit interaction on electronic 
states and it can be extended to study confined Ge struc- 
tures using an expanded basis of envelope functions^ 

This paper is organized as follows. Section [TXI discusses 
the symmetry of wavefunctions and their effects on spin- 
flip processes. Section Hill provides a quantitative discus- 
sion of intervalley spin flips due to scattering with zone- 
edge phonons. Using group theory, we derive the spin-flip 
matrix elements and calculate the resulting spin lifetime. 
A comparison with numerical calculations is also pro- 
vided. In Sec. IIVI the method of invariants^— is used 
to construct a spin-dependent k-p Hamiltonian around 
the L point. We solve this Hamiltonian to get the energy 
bands and spin-dependent wavefunctions. This informa- 
tion is then used in Sec. |V]to study intravalley spin flips 
due to scattering with long-wavelength acoustic phonons. 
Section [VII summarizes our findings. 
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FIG. 1: (a) The four conduction-band valleys of Ge. These 
valleys have ellipsoidal energy surfaces and their centers are 
located at the zone-edge L points, (b) Band structure of Ge 
along high-symmetry crystallographic axes. The conduction- 
band edge is indicated by the irreducible representation L\ 
(Le) in single (double) group notation . 49 ' 50 



II. SYMMETRY CONSIDERATIONS 

Spin flips during an electron scattering are described 
by the matrix element (k, -f|-|i/j|k', JJ.) where Hi denotes 
the scattering interaction and s = {f|-, .Jj.} are the spin-up 
and spin-down states. In a multivalley conduction band 
this scattering represents an intravalley or an interval- 
ley process depending on whether the electron wavevec- 
tors k and k' reside in the same or different valley (s), 
respectively. Figure [Ija) shows the four valleys in the 
conduction band edge of Ge. If the crystal possesses an 
inversion center then each band at wavevector k is spin 
degenerate and we can define its states with respect to 
the spin orientation n, 

(k,/i,ft |<r-n|k,/i,ft) = -{k,n,i\. \a-n\k,fiA) > , 
(k >l i,lMo-.n|k,/x,^)=0, (1) 

where ^ is the band index. In the following we omit this 
index and deal only with the lowest conduction band. 
Mixed by spin-orbit interaction, the spin-up and spin- 
down states (ft,J|) in Eq. ([1} are not pure spin states 
(t, I) • They are written as^ 1 - 



|k,1r> 



4Kt> + ^», 



E( a k*i^>-4>t>), 



(2) 



(3) 



where £ runs over a basis of periodic Bloch functions. Us- 
ing the crystal and time reversal symmetries at the four 
L-valley centers, we can identify a handful of periodic 
functions that describe the spin properties of conduction 
electrons in Ge. The analysis relies on their transforma- 
tion properties under symmetry operations. 

Representations of conduction (valence) states in the 
L point have even (odd) parity under space inversion op- 
eration. In the notation of single group theory, there 
are 6 irreducible representations (IR) of the L point 
space groupi 49 i 50 The lowest conduction band is non- 
degenerate and belongs to L\ (a one-dimensional IR). 
Figure [IJb) shows that nearby bands are pairs of va- 
lence and conduction bands. They are represented, re- 
spectively, by the two-dimensional IRs L' 3 and L 3 . If the 



crystal potential is vanishingly small then the L point en- 
ergies of these five bands (Li, L 3 and L' 3 ) are degenerate. 
However, the crystal potential in Ge splits these bands 
into three sets and the relatively large energy separation 
from L\ will be shown to result in a very slow intravalley 
spin relaxation process. In comparison, the six conduc- 
tion band valleys in Si are located close (in energy and 
wavevector) to the two-band degeneracy at the X point. 
This degeneracy leads to a spin hot-spot along certain 
directions in the square boundary faces of the Brillouin 
zone^ i 10 ' 13 and to a unique behavior of intravalley spin 
relaxation in Si. 13 This distinct difference between Si and 
Ge merits independent treatments of the spin relaxation. 

The wavefunctions of conduction electrons, |k,-|f (J|)), 
include small contributions from states of remote bands. 
In the kp theory there are two first-order terms related 
to signatures of the spin-orbit interaction, 



H SO 



h 



4toqC 2 



W x p ■ a, 



W x k • a. 



(4) 
(5) 



Hg transforms as a polar vector (W) and can couple 
between odd and even states. In our case, the coupling 
is between states of L\ and L' 3 symmetries that repre- 
sent, respectively, the lowest conduction band and upper 
valence bands. On the other hand, Hso transforms as 
an axial vector ( W x p) and can couple the even states 
of L\ and L3 (lowest and upper conduction bands). The 
dimensionality of L' 3 or L3 (2D IRs) is such that each is 
coupled to L\ by two components of a vector that lie per- 
pendicular to the valley axis [i.e., parallel to the hexago- 
nal boundary faces at the zone edge; see Fig.QJa)]. 

The aforementioned /c-independent coupling between 
Li and L3 zone-edge states can be used to find new eigen- 
states that diagonalize the L point Hamiltonian (k=0) in 
the presence of spin-orbit coupling, 



1^1. ->+ E E 

s'=t4 t 



iL^s'jiL^s'lHsolLus) 



Ee — El 1 
Li,s) +bx\L 3 ,s) + b 2 \L Zl -s) 



(6) 



where 61,2 -C 1 are fc-independent parameters, —s is the 
opposite spin direction of s, and \Li,s) denotes \Le}® |s). 
Slightly away from the valley center, the signature of L 3 
on conduction electronic states outweighs the signature of 
L' 3 via Hg (since k <C 27r/a where k is measured from 
the L point and a is the lattice constant). This prop- 
erty implies stronger spin-orbit coupling with the upper 
conduction bands than with the upper valence bands. 
Therefore, we start with a study of the intervalley spin 
relaxation where the effect of Hso enters at the lowest 
order (the spin-flip amplitude does not depend on rela- 
tive value of the wavevectors with respect to their valley 
centers). 
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III. INTERVALLEY SPIN-FLIP TRANSITIONS 

The intervalley scattering in Ge is governed by elec- 
tron interaction with zone-edge phonons close to the 
X point. For example, the L m and L n i valley cen- 
ters [k = 7r/a(l, 1, ±1)] are connected by the wavevec- 
tor Xqoi [k = 27r/a(0,0, 1)]. These zone-edge phonons 
belong to three 2D IRs: X 3 (TA, 10 meV), X x (LA & 
LO, 29 meV) and X 4 (TO, 33 meV)^i with their po- 
larizations and energies written in the parentheses. If 
the transition involves valleys that we denote by L and 
L t , then the selection rule for electron scattering reads 
Li <£) Lu — X\ © A^ where A3 is forbidden by time re- 
versal symmetry. 52 ' 53 

Including the spin degree of freedom, it is convenient 
to use double group theory where L\ is replaced by L§. 
The symmetry properties of the renormalizcd state in 
Eq. ([6]) are described by the L§ IR. The new selection 
rule readsi^ 



J 6(t) 



2A 1 0A 4 ©X., 



(7) 



and it means that three independent scattering param- 
eters are needed to fully describe intervalley scattering 
(two are related to the X\ symmetry and one to A4). 
Scattering with zone-edge phonons of A3 symmetry is 
excluded at the lowest order by time reversal symme- 
try. This selection rule does not provide information on 
the spin orientation dependence. To overcome this short- 
coming, we can work with IR matrices rather than their 
tracesji 3 - Technical details of applying this approach in 
Ge are given in Appendix A and below we present the 
main findings. 

Spin flips due to scattering with phonons of A4 sym- 
metry are described by a single independent nonvanish- 
ing matrix-element constant which we denote by Dx 4 , s - 
On the other hand, of the two independent nonvanishing 
matrix elements of X\ , one is attributed to spin- flip tran- 
sitions {Dx 1 . s ), and the other to spin-conserving scatter- 
ing {Dxx.m)- Most important in analyzing experiments 
is the dependence of the scattering amplitude on the spin 
orientation of electrons. The spin orientation, n, is de- 
scribed by a polar angle from the +z direction (9) and 
an azimuthal angle in the xy plane measured from the 
+x direction (</>). For a spin-flip transition between Lm 
and L n y valleys (via a zone-edge phonon with wavevec- 
tor X001), the square amplitude of the matrix element 
reads 



|(k Lt ,4„|i^|kL^„)| 2 

J=l,2 

2D Xi s {l + cos 2 9 + sin 20 sin 2 1 



(8) 



sin 



if % 
if % 



j sums the two degenerate branches of each X point IR. 
For the remaining five transitions between other combi- 
nations of L points, Eq. ([8]) varies according to a straight- 
forward coordinate transformation. These results are 
summarized in Table Q] (also see discussion of Eq. (|A23[) 



TABLE I: E j=a ,b \ M xl ( k ^> k "' ^)/2D x , | 2 for interval- 
ley spin-flip transitions of all the six valley-to-valley configu- 
rations. 



L L t 


A'i 




X 4 


tin <-> L L1 - 


1 + cos 2 6 + sin 2 8 sin 2(f> 




1 — cos 8 


Liu L 1 j 1 


1 + sin 2 8 sin 2 <f> + sin 28 cos (f> 


1 


— sin 2 8 sin 2 (f> 


Lm Ljjj 


1 + sin 2 8 cos 2 <f> + sin 28 sin (f> 


1 


— sin 2 8 cos 2 <j> 


Lyu <-> L llT 


1 + sin 2 8 sin 2 <f> — sin 28 cos 4> 


1 


— sin 2 8 sin 2 4> 


Lj-11 -^lTl 


1 + cos 2 8 - sin 2 8 sin 2</> 




1 - cos 2 8 


Ljjj L n j 


1 + sin 2 8 cos 2 4> — sin 28 sin (f> 


1 


— sin 2 8 cos 2 (j> 



in Appendix A). As seen by the angular dependence in 
Eq. ([5]), spin-flip transitions are forbidden for specific 
choices of the spin orientations. For example, if the spin 
is oriented parallel to the phonon wavevector (z-axis in 
this case and 9 — 0), then the spin-flip amplitude van- 
ishes for the A4 symmetry. Below, we will quantify the 
effect of spin orientation on the spin lifetime. 

Having the spin-flip matrix elements, one can readily 
calculate the intervalley spin relaxation rate by integra- 
tion over the Brillouin zone, 



1 



2nh 

Wc 



d A k 



df(E) 



dE 



d 3 k' |M^(k,k')| 
(2^)3 n„(q) 



(9) 



2(n„, q + i ± \)5{Eu -E k ± n v (q)). 



q = 5.323 g/cm 3 is the crystal density, f{E) is 
the statistical distribution of electronic states, and 
N c = J d 3 kdf(E)/dE\ E is an effective density of states 
constant, v, q = k k', Sl^q) and n^.q denote, re- 
spectively, the phonon mode, its wavevector, energy and 
Bose-Einstein distribution. The ± refers to phonon emis- 
sion and absorption processes. M^(k, k') = (k',J| 
|M,/ q|k, ft) is the spin-flip matrix element. 

Using the spin- flip matrix elements in Table U we cal- 
culate the intervalley spin relaxation rate for a Boltz- 
mann distribution of electrons, 



1 



T~s inter 



2md>? 



E 

;=i,4 



A, 



h 2 Q^/TTi exp(y l )-l 



(10) 



Tdd— y mirrit is the effective electron mass where 
rat ~ 0.08mo and mi rj 1.6mo are, respectively, the trans- 
verse and longitudinal components. r &(yi=fti/kBT) = 
y / y7exp(?/i/2) K-i(yi/2) is associated with the modified 
Bessel function of the second kind. This term slightly 
depends on temperature and varies between 2 and 4 in 
the temperature range between 10 K and 400 K (for 
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both phonon energies). On the other hand, most of 
the temperature dependence of the intervalley relaxation 
rate comes from the thermal population of zone-edge 
phonons (exponent term in the denominator) . This pop- 
ulation is strongly suppressed at low temperatures. Fi- 
nally, the scattering constants are Dx 1 , s = 35meV/A 
and Dx 4 ,s = 46meV/A, quantified by numerical results 
of empirical pseudopotential method; 54 adiabatic bond- 
charge model^ 5 - and rigid-ion approximation^ They are 
respectively weighted by Ai(0,4>) and Aj(0, </>) that in- 
clude the dependence on the spin orientation. We discuss 
their values for several general cases. 

no-strain or [100] -strain: The four L- valleys in the 
lowest conduction band are degenerate and transitions 
between all six pairs of valleys are equivalent. The 
anisotropy in spin relaxation due to intervalley scat- 
tering between two valleys is compensated by opposite 
anisotropy of other pairs. The sum of expressions in each 
of the two columns of Table fl] is independent of 9 and </>, 

Ai=8, At =4. (11) 

As shown next, when the symmetry between different val- 
leys is broken, the dependence of the intervalley matrix 
elements on spin orientation lends itself to a measurable 
anisotropy in the spin lifetime. 

[Ill] strain: The case of uniaxial compressive strain 
results in a single low-energy valley (along the strain axis) 
and three higher energy valleys. At relatively large strain 
levels (~1%), the energy split is large enough to quench 
the intervalley spin relaxation mechanism^ This effect 
amounts to assigning 1/T s> i n t er — 0. On the other hand, 
in biaxial compressive strain configuration (or uniaxial 
tensile strain) three of the valleys shift down in energy 
and one valley shifts up. Excluding transitions with the 
Lux valley (considering the last three lines in Table QJ 
we get 

16-4sin 2 6'sin2(/>-4sin26'(cos(/)+sin0) 8 
A l = jj , ^4=g - (12) 

This strain configuration restores the anisotropy in spin 
relaxation due to electron scattering with X\ zone-edge 
phonons. By changing the spin orientation from the 
strain axis to its perpendicular plane, r Sj i nter drops by 
-50% [changing A x from 8/3 to 20/3 in Eq. JTD])] . 

[Hi)]- strain: This strain configuration is optimal for 
detection of the anisotropy since intervalley transitions 
are effective from a single pair of valleys. Consider, for 
example, the case that Lm and L n y valleys shift suffi- 
ciently down in their energy. Then, only the first term in 
Table U represents the intervalley scattering and we get 

Ai = 2(l + cos 2 6> + sin20sin 2 6>), A 4 = 2sin 2 6». (13) 

The anisotropy in spin relaxation is now caused by elec- 
tron scattering with both types of zone-edge phonons. 
By changing the spin orientation from the strain axis to 
its perpendicular plane, T s ^ nter [Eq. (|10p ] is doubled. 




10 50 100 200 10 100 200 300 400 
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FIG. 2: Temperature dependence of the intrinsic spin life- 
time in unstrained Ge for spins oriented along the z- 
crystallographic axis. The solid lines follow Eq. (|10[) and 
denote the spin lifetime from intervalley scattering with X\ 
and X± phonons [the temperature range is extended in (b)]. 
At room temperature, the spin lifetime of this scattering is 
~ 1 ns. The markers in (a) are numerical results. Black di- 
amonds denote contributions from X3 phonons and red pen- 
tagrams from X\ and X4 phonons. The green circle mark- 
ers in (b) show the combined effect. The dash red line in 
(b) denotes the spin lifetime from intravalley scattering with 
long- wavelength acoustic phonons [Eq. (|33|) ]. Yellow square 
markers are the respective numerical results. 



Other than strain, it should also be possible to observe 
the anisotropy by applying electric fields of few kV/cm 
along the mentioned directions. Here, valley repopula- 
tion will result in preferential scattering from hot-to-cold 
valleys j££ Finally, by averaging over spin orientations 
in Eqs. ([12 ]) - ([T5 ]) . the spin lifetime [Eq. ([10])] with two 
(three) low-energy valleys is about 3 (3/2) times longer 
than that of the unstrained case. The reason is that elec- 
trons can scatter to one (two) valleys rather than three. 

The solid blue curve in Fig. HJa) shows the temper- 
ature dependence of the intervalley spin lifetime in un- 
strained bulk Ge [Eq. (USD with Ai=8 and A 4 =4]. We 
have also performed numerical integrations of Eq. Q 
following the procedure in Ref. Q. The numerical cal- 
culation of the electron-phonon matrix element follows 
a rigid-ion approximation! 56 ' 58 Whereas this numerical 
approach is not transparent compared with the group 
theory analysis, it takes into account higher-order cor- 
rections (variation of the matrix element when departing 
from the center of the valley). Nonetheless, the numerical 
results for scattering with X\ and X 4 phonons [red pen- 
tagram markers in Fig. H£a)] show that the zeroth-order 
analytical calculation is sufficient [Eq. (ITUl) ]. The black 
diamond markers denote numerical results due to scatter- 
ing with X3 phonons. While their zeroth-order contribu- 
tion vanishes by time reversal symmetry [Eq. ([?])], their 
first-order contribution (linear in wavevector) becomes 
important at low temperatures. The reason is that their 
low energy (—10 meV) leads to a much larger population 
at low temperatures compared with that of X\ and X4 
phonons (whose energies are ~30 meV). 
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IV. L POINT HAMILTONIAN 

In this section, we expand the wavefunctions of elec- 
trons using the L point basis states. Important signa- 
tures of spin-orbit coupling on these wavefunctions will 
be identified and correlated with intravalley and interval- 
ley spin relaxation processes. Near the valley center, the 
wavefunction is approximated by 



Ik,s) = 



[ £ 

7 = 1,3,3' 



C 7 (k, s ) 



i(ki+k)-r 



(14) 



where following Eqs. ([S])-©, we write 



JV-, 



C 7 (k,ft)|L 7 ) =J2 a m , 7 (k)|L7, t) + b m ^k) \L™,1), 

771= 1 

N-, 

C 7 (k,4)|L 7 > =J2<,4V\L™,1) - 6^ i7 (k)|L™,t>. (15) 



Totally, we consider 10 spin-dependent basis states: two 
from the non-degenerate lowest conduction band (N±=l), 
and four from either the upper conduction or valence 
bands (each being two-band degenerate in the absence 
of spin-orbit coupling, N 3 =N 3 >=2). The coefficients are 
eigenvectors of the Hamiltonian, 



H 33 + E, 




d \=E Cj I, (16) 
C3/ \ C 3 / 



where Hij is a matrix block denoting the spin and 
wavevector dependent coupling between basis states with 
Li and Lj symmetries. E g ^ u and E g>v denote, respec- 
tively, the L point energy separations of the lowest con- 
duction band from the upper conduction and upper va- 
lence bands [see Fig.QJb)]. Below, we present the Hamil- 
tonian matrix using the basis functions of the Lni point 
[k L _= tt(1, 1, l)/a]. Matrix forms in the (111), (111), and 
(111) valleys are derived by trivial coordinate transfor- 
mation. In addition, to derive a compact matrix form we 
use a rotated set of Cartesian coordinates, 



x — y _ x + y — 2z „ x + y + z 

W = — , U = — , A = 



V2 



V6 



v/3 



(17) 



u and w lie parallel to the hexagonal boundary face 
[Fig. QJa)]. a is along the valley axis connecting the T 
and Lin points. 

We construct the Hamiltonian matrix [Eq. (fT6|) ] using 
the method of invariants^ - — Application of this method 
with relevance to the L point is given in Appendix B. 
Here we summarize the findings. The lowest conduction 
band is associated with the identity IR and contributes 
a trivial 2x2 matrix form (L\ <g> L\ = Li), 



H n = 



2m* t 



2m* 



12x2 ■ 



(18) 



m* t and m\ are effective mass parameters representing 
the effect of remote bands (outside the chosen basis 
states). Matrix blocks of the upper valence bands or 
upper conduction bands share a similar form (L 3 ® L 3 = 
L' 3 ®L' 3 =L X +L 2 + L 3 ), 



h 2 kl 



2m* 



2m 



)/ 4 x4 + Aip y (g>a A , (19) 



where i — - 3 or i = 3'. The mass parameters have similar 
meaning as in Hn. Aj denotes the internal spin-orbit 
coupling between the two Li basis functions.— p y = <j v 
originates from the two-band degeneracy in the absence 
of spin-orbit coupling. 

The off-diagonal matrix block f/13/ denotes the cou- 
pling between the lowest conduction band and upper va- 
lence bands. Its form follows from L\ ® L' 3 = L' 3 , 

H 13 , = p(k w [0,l}-k u [l,0])®I 2H2 (20) 
+ ia (k x a) w ® [0, 1] - (k x <r)„ ® [1,0] , 

where [1,0] and [0, 1] are ordinary 1x2 matrices. Their 
Kronecker products with 2x2 matrices indicate that 
Hi 3 > is a 2 x 4 matrix. P and a are two independent ma- 
trix element constants that originate from the k • p and 
Hg perturbation terms, respectively. The coupling ma- 
trix of the lowest and upper conduction bands is wavevec- 
tor independent and it follows from L\ ® L 3 = L 3 , 



H13 



iA 



,(ov»® [1,0] 



[0,1] 



(21) 



where A^ denotes the direct spin-orbit coupling between 
these bands. Finally, the 4x4 coupling matrix between 
the upper valence and conduction bands follows from 
L' 3 ®L 3 =L' 1 +L' 2 + L' 3 , 

H 33 ' = [Pi(ik u p v - k w I 2x2 ) + P 2 k A p x ] ® I 2x2 , (22) 

where we have neglected the Hg coupling between these 
bands since it plays a negligible role in the spin relaxation 
of conduction- valley electrons. Table [V] in Appendix B 
lists the values of all parameters in Eqs. (IT8l - (|22l) . 

Given the relatively large L point energy gaps, the 
energy dispersion of electrons in L-valleys is well ap- 
proximated by eigenvalues of Hn + Hi 3 iH{ 3 , / E g . v — 
Hi 3 H\ 3 /E g , u , 



2A| 



h 2 (k 2 u + ki) h 2 k 2 



2m t 



2mi 



(23) 



The constant energy shift is due to the direct spin-orbit 
coupling with the upper conduction bands. The effective 
mass parameters are 



1 



1 

m * f 



2P 2 + 2a 2 
h 2 E„ v 



1 

mi 



1 

ml 



4a 2 



H 2 E g , v 



About half of the anisotropy between the transverse and 
longitudinal effective masses in Ge (m t w 0.08mo and 
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to/ ss 1.6too) is set by the spin independent coupling with 
the upper valence bands (P = 9eV • A). The spin-orbit 
coupling signatures on the energy dispersion are negligi- 
ble and can be ignored (a=40 meV-A and Al=27 meV). 
On the other hand, the minute effect of spin-orbit cou- 
pling on the eigenvectors of Eq. f|16 j) sets the timescale 
for spin relaxation. Choosing the spin quantization along 
the valley axis, the spin-up eigenvector along this direc- 
tion [n = A in Eq. (JTJ] reads 



Ci(k,ft A 
C 3 (k,ft A 

C|(k,ft A ) = 



0{k 2 ), 



E, 



g.u 
P 



0, -1, 0, i 



0(k 2 



E„ 



k u - ilskw, /+(k), 
k w - 173 k u , /-(k) 



0{k 3 ). (24) 



The components of the spin-down eigenvector [Cj(k, JJ. A )] 
are readily obtained from space inversion and time rever- 
sal relations [Eq. (|15|)]. The g(k) and /±(k) functions in 
Ci and C3 read 



?(k) 



P 2 

g,v 



/±(k) = r± [j!(k w - ik u ) ± i^ 2 k A ] 



(25) 



where r + = 1 and r_ = — i. The jj <C 1 parameters scale 
with three of the spin-orbit coupling constants (a, A3' 
and Al) 



7i = 



72 



73 



A L 


Pi 


Eg^u 


P 


a 


Al 


P + 


Eg :U 


a 




P + 


Eg t v 



0.006, 



ft 



0.005, 



0.05, 



(26a) 
(26b) 
(26c) 



where the internal spin-orbit coupling in the valence band 
(A 3 ') sets most of the value of 73. Only when the spin is 
oriented along the valley axis (n = A) , this parameter is 
excluded from the opposite-spin components of C3 (k, ft n 
) [i.e, from the /±(k) terms in Eq. (|24p]. This behavior 
will have important consequences on the anisotropy of 
intravalley spin relaxation. 



Connection between the L point Hamiltonian 
parameters and spin relaxation 

To facilitate a connection between the Hamiltonian 
eigenvectors and spin relaxation we make use of spin- 
flip overlap integrals. We show that the direct spin-orbit 
coupling between the conduction bands (Al) plays a key 
role in setting the intervalley spin relaxation rate (inde- 
pendently derived in Sec. lIII[) . On the other hand, we will 



see that intravalley spin- flip transitions are weaker. To 
make these connection clear, we write the overlap integral 



X(kU n ; k,ft n ) 



|Ct(k',4n)C 7 (k,ft n )|L 7 , kl ), (27) 



where k and k' are measured from the nearby val- 
ley center (k^ and k L >). The bra and ket states of 
this overlap integral include only the periodic Bloch 
parts in Eq. (fT4|) . While the combined phase factor, 
exp{i(k L — k L ' + k — k') • r}, is excluded from the over- 
lap integral, it will be taken into account in the phonon 
phase when calculating the matrix elements. Using 
Eq. ([M]) . the overlap integrals of electrons in different 
valleys reads (k L ^ k £ ') 

2" e (k',tr„ ; k,tr„) = cx. n 

I e (k',^ n ;k,^ n )=c 3 , n -^+0(fc 2 ), (28) 

where Cj n are constants of order unity which denote 
contributions from the spin-orientation dependence (n) 
and from the overlap of conduction basis states in dif- 
ferent valleys: (L± t k L \Lj^ L ,). Eq. (|2"8")) implies that 
the ratio between spin and momentum relaxation rates 
due to intervalley scattering is about A L /E 2 U (indepen- 
dent of the values of the wavevectors with respect to 
the valley centers). For intravalley scattering (k i =k £ ), 
on the other hand, the basis functions are orthogonal: 
(-^V.kJ-^k^/) = S u -fS mn . As a result, the spin-flip over- 
lap integral for electrons of the same valley reads 



2P 2 
E 2 

9,v 



{"/2q+K A + i"fiq-K-), (29) 



where q = k-k', 2K = k+k', and X± = X w ±iX u . The 
terms have quadratic wavevector dependence and they 
are proportional to the spin-orbit constants in Eq. (|26|). 60 
The overlap integral of other spin orientations (n ^ a) will 
be discussed in the next section. 



V. INTRAVELLEY SPIN RELAXATION 

The power-law dependence of intravalley spin-flip ma- 
trix elements can be identified by their transformation 
properties under time reversal and space inversion op- 
erations. Yafet showed that spin-flip matrix elements 
due to scattering with long- wavelength acoustic phonons 
have a cubic (quadratic) wavevector dependence in Ge 
(Si)4i In Appendix C, these important findings are gen- 
eralized and it is shown that in Ge (k',JJ- |iJj^ tra |k, ft) 
scales with Kiq m q n for scattering with acoustic phonon 
modes (A = LA or TA) and with Keq rn for optical phonon 
modes (A = LO or TO). K and q are, respectively, the 
average and difference of k and k'. For intravalley scat- 
tering in Si the K dependence drops. We first explain 
this interesting difference. 
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From inspection of the wavevector dependence of in- 
travalley spin flips in Ge (Kq l ), one sees that they are 
forbidden between opposite points with respect to the 
valley center (K = 0). This restriction on spin- flip tran- 
sitions is a manifestation of time reversal symmetry. In 
silicon, K-dependent scattering belongs to the intervalley 
g-process which involves transitions between two valleys 
on opposite sides of the same crystal axis i 10 i 13 Since in 
Ge the valley center is at the zone edge (L point), this 
type of scattering occurs within a single valley. Its de- 
pendence on the wavevector components (Kq l ) amounts 
to the combined effects of intervalley g-process and in- 
travalley scattering in Si (K and q l ). 

Beyond the power-law dependence, an analytical ap- 
proach to derive accurate intravalley matrix elements 
requires a combination of k ■ p, rigid-ion and group 
theories.— Because of the wavevector dependence of 
these matrix elements, one cannot invoke group theory 
alone to find their exact forms (as we did for zeroth-order 
intervalley spin flips). We employ a simpler approach 
than in Ref. [l3[ and describe the interaction with long- 
wavelength acoustic phonons by H™/^ A = Eq where 3 
is an effective deformation potential constant ! 10 ' 49 This 
scalar form averages out the scattering angle dependence 
of the second-rank deformation potential tensor 4^ We do 
not model the electron scattering with long-wavelength 
optical phonons since it is a weak effect in nonpolar 
semiconductors j 49 ' 61 

We use selection rules of the L point space group to 
construct the spin-flip matrix element from the overlap 
integral. The transformation property of the deformation 
potential tensor, L3 (8-E/3 = L1+L2 + L3, implies that di- 
rect coupling of conduction and valence states is excluded 
because of their opposite parities {L\ ® I/3 = L'3). This 
tensor can, however, couple any of the basis states to 
themselves (Lj <8> Li). This behavior justifies the use of 
the spin-flip overlap integral. The resulting intravalley 
spin-flip matrix element in the Lm valley is approxi- 
mated by 



M°f tra (k>, k; n) a Zql a (k> , ^ ; k, ft n ) • 
Following a straightforward procedure we find 



(30) 



X (k',JJ. n ; k,fr„) =isintfA t +A l cos 2 - + A\ sin 2 - ,(31) 
where cos # = n • A and 



1 

2' 
2P 2 



1, = -(Z a (k',^ A ; k,tO-X a (k',^ A ; k,^ A )) 
-73 (K X q) A , 



E 2 

g,v 

Ai = e -^i a {k',U; KU) 

2P 2 



E 2 

g,v 



(72 g+ if a + ijiq-K-)t 



(32) 



<p is the azimuthal angle of n measured with respect to 
the w-axis in the u>u-plane. Most importantly, 73 which 



incorporates the effect of the internal spin-orbit coupling 
in the valence band [Eq. (|26cj) ] does not affect the spin- 
flip amplitude [Eq. (I30[) ] when the spin orientation is 
along the valley axis (1? = 0). This effect will lead to 
a pronounced anisotropy in the intravalley spin lifetime. 

It is not surprising that the overlap integral approach 
yields correct wavevector power-law dependence [substi- 
tuting Eqs. ([3"I ]) -(|3"2" ]) into Eq. flUJ)]. The space inver- 
sion and time reversal symmetries are respected by the 
Hamiltonian whose eigenvectors were used to calculate 
the intravalley spin-flip overlap integral. These symme- 
tries also lead to the so-called Elliott- Yafet cancelation of 
all terms up to quadratic order in qj^dl. In fact, since the 
Hamiltonian respects all other symmetries of the L point 
space group, the intravalley matrix element shows other 
selection rules^ 2 - From Eq. (|3"2"|) we see, for example, that 
a spin-flip is forbidden when the electron is scattered 
along the valley axis (i.e., q w =q u =0, q A y^0). This con- 
straint is understood by the symmetry of the vector-type 
coupling with the valence states {Li®L'% = L'3). As men- 
tioned, this coupling involves the two transverse compo- 
nents (w and u) with respect to the valley axis (a). 

We calculate the spin lifetime in the Lm valley due to 
electron scattering with long- wavelength acoustic phonon 
modes. This intravalley process dominates the spin re- 
laxation under conditions of [111] strain. 1 ^ For sufficient 
uniaxial compressive strain along this direction (~ 1%), 
one of the valleys is significantly lowered in energy and 
the intervalley scattering is quenched. Then, phonon- 
induced intravalley spin-flips can dictate the spin relax- 
ation of conduction electrons if scattering from impuri- 
ties is negligible (non-degenerate doping). To get an an- 
alytical expression of the intravalley spin lifetime, the 
phonon energy is approximated by ^Ac(q) — fov AC q 
where v AC — 3.5 • 10 5 cm/sec is the speed of acoustic 
phonons in Ge. We also make use of the long- wavelength 
limit and approximate the acoustic phonon population 
by fcsT/f2 AC (q) 3> 1. Then by considering a Boltzman 
distribution of electrons and substituting Eqs. (|30l) - (|3"2")l 
into Eq. ([9]) one gets 



1 _ll 

'8, intra Tq 



k B T 
U 



[sin 2 + (1 + cos 2 , (33) 



where Uq = 25.8 meV is the room-temperature thermal 
energy, « 0.12 and r « 0.3 ns are expressed by 



P 



2mal + 3m t j 2 



hm t % 



1 



1024 



E 2 \2t:J h 4 gv 



2 ' 

AC 



(34) 
(35) 



In accord with momentum scattering, we have used 
a value of S = 7.5 eV for the deformation potential 
constant 4^ The anisotropy in the intravalley spin relax- 
ation is evident [square bracket term in Eq. (|33[) ]. Our 
analysis shows that the lifetime is the longest for spin ori- 
entation along the valley (and strain) axis where $ = 0. 
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It drops by nearly a factor of 5 when the spin is oriented 
in the perpendicular plane (# = 7r/2). At room temper- 
ature, this change amounts to reducing the intravalley 
spin lifetime from ^700 ns to ^150 ns. These extremely 
long timescales are a consequence of the space inversion 
symmetry and the position of the valley center in the 
edge of the Brillouin zone. 

The temperature dependence of the intravalley spin 
lifetime is shown by the dash line in Fig. Hfb) for spin 
orientation along the z-axis [assigning cos 2 $ = 1/3 in 
Eq. (|33[) ]. In this spin orientation, the intravalley spin 
relaxation rate is equivalent in all four valleys. Fig- 
ure [2Kb) also shows that in unstrained bulk Ge, the spin 
lifetime of conduction electrons due to intravalley scat- 
tering with acoustic phonons is two orders of magnitude 
longer than the intervalley spin lifetime at room temper- 
ature. In addition, at very low temperatures the intrinsic 
spin lifetime reaches timescales of one second. Therefore, 
the phonon-induced spin relaxation is likely to be read- 
ily masked at low temperatures by localization effects on 
residual impurities (e.g., hyperfine interactions and Ra- 
man processes) ji 3 - Finally, the square markers in Fig.[2{b) 
show results of rigorous numerical calculations following 
the procedure in Ref. Evidently, the overlap integral 
analytical approach provides rather accurate results and 
yet it clearly explains the underlying physics. 

Before concluding this part, we compare three aspects 
of the intravalley spin relaxation in Si and Ge. First, the 
overlap integral approach is valid in Ge due to the rela- 
tively large separation of the non-degenerate conduction 
band from other valence and conduction bands. In Si, 
on the other hand, the intravalley spin relaxation is af- 
fected by the proximity of the conduction bands where 
the off-diagonal terms of the deformation potential play a 
key rolei 10 ' 13 Second, along the A-symmetry axis which is 
relevant in Si, the spin-orbit coupling does not lift the en- 
ergy degeneracy between the upper pair of valence bands. 
As a result, the intravalley spin relaxation is not affected 
by the internal spin-orbit coupling in the valence band 
and the anisotropy is weaker in Si reaching a factor of 
twoJ^ Finally, the intravalley spin relaxation rate ex- 
ceeds the intervalley rate below 50 K in SijS 3 - and below 
20 K in Ge. Reasons for the difference are the larger en- 
ergy of zone-edge phonons in Si and the T 5 ' 2 rather than 
T 7 / 2 temperature dependence of its intravalley process. 



VI. SUMMARY 

We have presented various origins that limit the in- 
trinsic spin lifetime of conduction electrons in Ge. In 
unstrained bulk Ge and at T > 20 K, the intrinsic spin 
lifetime is limited by intervalley electron scattering with 
zone-edge phonon modes of X\ and symmetries 
(reaching ^1 ns at 300 K). This spin lifetime is gov- 
erned by the coupling with the upper conduction bands 
and its temperature dependence is set by the thermal 
population of the zone-edge phonons (with energies of 



about 30 meV) . By analyzing the crystal and time rever- 
sal symmetries in the multivalley conduction band, we 
have found the spin orientation dependence of the domi- 
nant intervalley spin-flip processes. This dependence al- 
lowed us to quantify the change in the intervalley spin 
lifetime when varying the spin orientation under various 
stress configurations [Eqs. (jTU| - (fT5|) ]. 

We have derived a spin-dependent kp Hamiltonian 
model in the vicinity of the zone-edge L point [Eqs. (|16[1 - 
(|2"2"1) ]. This compact model provides a lucid picture of 
the spin-orbit coupling effects in Ge. Similar to using the 
Kane model in zinc-blend semiconductors^ the compact 
L point Hamiltonian has implications beyond derivation 
of spin-flip matrix elements. For example, by employ- 
ing a plane-wave expansion along confined directions in 
nanostructures, this Hamiltonian model can be used to 
study spin-dependent properties in Ge nanostructures. 
Together with Si related theories ; 10 ' 13 one can also inves- 
tigate spin properties in SiGe alloys. 

Using the eigenvectors of the Hamiltonian matrix, we 
have derived forms of the spin-flip matrix elements due 
to intravalley scattering with long-wavelength acoustic 
phonons [Eqs. (|30l) - (|3"2"|) ]. The intravalley spin lifetime is 
found to be two order of magnitude longer than the in- 
tervalley spin lifetime. As such, intravalley spin-flips af- 
fect the overall spin relaxation only when quenching the 
intervalley spin relaxation (e.g., by application of a uni- 
axial compressive stress along the [111] crystallographic 
axis)*^ 2 - Beyond the T 7 / 2 temperature dependence of the 
intravalley spin relaxation, we have also quantified its 
dependence on the spin orientation [Eq. ([33)1]. The 
anisotropy of the intravalley spin relaxation results in 
a remarkably long spin lifetime (nearly 1 (is at room- 
temperature) when the spin is oriented along the val- 
ley (and strain) axis. The relatively large anisotropy 
of the intravalley spin relaxation was explained by the 
coupling with the internal spin-orbit interaction in the 
valence band. 

We have elucidated the differences in the spin relax- 
ation of bulk Si and Ge crystals. While both materials 
have a diamond-crystal structure, in Ge the valley center 
is located at the edge of the Brillouin zone (L point) and 
the lowest conduction band is well separated from other 
bands. These properties lead to a very long intravalley 
spin lifetime in Ge with a cubic power-law dependence 
of intravalley spin-flips on wavevector components. This 
cubic dependence is also expected to be larger than in 
graphene where unlike Ge but similar to Si, the time- 
reversal operation couples states in inequivalent valleys. 
Therefore, in spite of being heavier than Si and carbon, 
non-degenerate and strained bulk Ge is a very promising 
material choice for implementing spintronic devices^ - — 

This work is supported by AFOSR Contract No. 
FA9550-09-1-0493 and by NSF Contract No. ECCS- 
0824075. 



9 



Appendix A: Derivation of the selection rules for 
intervalley spin-flip transition 

We first focus on scattering between the Lm and L u j 
valley centers [k L = (1, 1, l)/2 and k Lt = (1, 1, -l)/2]. 
Generalization to other valley centers is made at the end 
of the Appendix. 

The selection rules connecting L and L t points involve 
common symmetry operations of the little groups at k^, 
-k Lt and q x = k Lt - k L , 



TABLE II: Non-trivial relevant IR characters and matrices in 
a intervalley scattering between kL and k_L t valleys. For ID 
IR ii(t) an d 2D IR X3, only characters are used and shown. 
/ and a x used in 2D IR X\ and X4 are the 2x2 identity ma- 
trix and Pauli matrix. These matrices are based our choice 
of basis states. The final results do not depend on this spe- 
cific choice since the two phonon modes belonging to each IR 
are degenerate. xZ t Lt — Xl" • Also shown is the effect of 
exchange operations on L star. Basis states in D 1 / 2 is along 
±z in spin space. 



g c E { (e\0),(e\0),(32 X y\r),(S 2x y\r), 

(i\r),(i\r),(p xy \0),(p xjj \0) }. (Al) 

They also involve operations that switch between k^ and 
-kit, 

g e e { (8 2z \0),(S 2z \0),(Pz\t),(Pz\t), 

(p xy \0),(p X y\0)},(6 2xy \T),(62 X y\T) }. (A2) 

The bar over operations denotes an additional 2tt rota- 
tion (in double group notation). Table HTl lists the char- 
acters of the nontrivial operations. By considering these 
operations and time reversal symmetry, the number of in- 
dependent nonvanishing matrix elements for each of the 
zone-edge phonon symmetries in diamond-crystal struc- 
tures (Xi, A 3 , A 4 ) is given by 



A/'; 



A", 



1 

2h 



2>£ 



t (5c)x^ + (ffc)x^(5c) 



(A3) 



where ho = 8 is the number of g c or g e operations and 
X L + = XL 1(t) x Xi/2- The second sum in Eq. (|A3|) de- 
notes the effect of time reversal symmetry and the minus 
sign takes into account the parity from the spinor basis 
and interaction H ep (see, Ref. |45[ for more details). By 
straightforwardly plugging the characters of Table HI1 into 
Eq. (|A3I) one finds the general selection rule of Eq. J7]) . 

Our aim is to express interaction matrix elements 
(kLt, S2\Hxi |k_L, si) between specific spin species in 
terms of Af Xi independent constants. In order to identify 
these scattering constants, we connect different matrix 
elements via appropriate symmetry operations. 

First, by time reversal and space inversion symmetries 
of diamond-crystal structures we can write 

(k Lt A\H Xl \k L ,H) = (k Lt A\H Xl \k L Ar, (A4) 
(k Lt A\H Xl \k L ,H) = -(k Lt ,^\H Xi \k L Ar-(A5) 

These identities hold for all phonons and possible spin 
orientations. The minus sign in Eq. (IA5|) roots from the 
Pauli matrix a y in the time reversal operator T = Ka y , 
where K is the complex conjugate operator. 



(S 2 xy\r) 



(i|r) 



(fe|0) 



{Pz\t) 



(Pxy\0) 



{&2 X y\r) 



X, 



-I 



X, 



Li 



Lu 



kiO-kit 



ki,o-ki, t 



kto-kit 



kiO-kit 



/2 



e . ¥ /0 , 



1 cf 

1 



-i N 
i 



-i Cf 
i 



<n 1 



We first study the case n||z, where n is the spin ori- 
entation. For Ai, the (p X y\0) operation equates spin- 
conserving transition to itself (seen from the IR matrices 
ofAiand£> r + = D Lut , x D in in Table©, 

6(t) > 1 

(k it , H- \H Kib) \k L A) (P = 0) (ki t ,^ |fl Jc . w |k i ,4),(A6) 

where two Ai basis states are denoted as Af and X\. 
This choice is arbitrary and will not affect the final re- 
sults due to the 2-fold degeneracy of the phonon modes. 
Eqs. (|A4|) and (|A6[) require the matrix elements of each of 
the Ai phonon branches to be a real number. From Ta- 
ble mi one can also find (i\r) relates the matrix elements 
of the two degenerate modes by a minus sign, 



(k^ ft l-ffy.wlki.1r) 



(i\r) 



(kit,ft|ff xi , w |ki,ft).(A7) 



With this additional information, a real number D Xlim 
could be assigned such that 

(kit, ft|ff* f |ki, ft) = -(k Lt , ft|ff xf |ki, ft> =D Xl>m ,(A8) 

Other operations do not give further information on these 
matrix elements. 

With the same operations, the result for spin-flip tran- 
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TABLE III: £ \M x j(k L ,ft n ;-k Lt An)/D Xi \ 2 for intervalley spin flips between Lm and L n j valleys. For each of the non- 
vanishing modes, Xi, the relative amplitude is provided for spin orientation (n) along any of the inequivalent high-symmetry 
crystal directions. Results between other valleys can all be obtained by trivial rotation transformation. 



n 


[0 1] 


[10 0] 


[110] 


[110] 


[10 1] 


[111] 


[Hi] 


Xi 


4 


2 


4 





3 


4 


4/3 


x 4 





2 


2 


2 


1 


4/3 


4/3 



sition is 

(k Lt A\H xam \k L , 1tf = 0) -i(kLu it\H x a W \lL L , (A9) 



(k Lt A\H xa{b) \k L A) 



(i\r) 



k it ,^|if ;c0(i , ) |k i) f).(AlQ) 



Together with Eq. (|A5|) . we can assign a real number 
Dx x ,s such that 

(k Lt ^|i^ f |k L ,^) = -(k Lt A\H xbi \k L A) 

= (l + i)D XuS . (All) 

Next we analyze matrix elements due to modes, 
where there is only one independent scattering constant. 
From Table HH the operations (52xy\T) and (i\r) give re- 
lations for spin-conserving transitions 

(k Lt ,it\H <b) \k L Af 2 = T) -(k L tA\H x:W |k L , ^),(A12) 
(k Lt ,1t\H xam \k L ,1t) ( = ] -<k it ^|ff-., (a) |k i ^).(A13) 



Together with Eq. (|A4|l . a real number D Xi . s could be 
assigned 

(k Lt ,it\H Xi \k L ,it) = -(k Lt ,it\ H xl \ k L ,H)=iD Xi , s .(AU) 

For spin-flip transitions, the exchange operation (c>2z|0) 
together with the general time reversal operation connect 
the matrix elements to their negatives, 



k Lt A\H xa{b) \k L ,H) 



(A15) 



(-k L A\H x:ib) \-k Lt A) T J -(k Lt A\H x: < b) \k L A), 

where the time reversal operation sends electron states to 
their Kramers conjugate, and keeps the electron-phonon 
interaction. Thus spin-flip matrix elements duo to Xi 
phonon modes vanish. 

Therefore, with spin direction along z, the scattering 
matrices from kL to ku for relevant phonon modes are 



, D xi,m {-l + i)D Xlt 
Hx? = -H xb =\ |,(A16) 
(1 + i)D Xus D x m 



H x? = ~H X b = 



iD 



X 4 ,s 







-iD Xl>s , 



(A17) 



where Eqs. (|A4j) and (|A5j) are used to get two other el- 
ements in each matrix. Eqs. (|A16j) and ()A17|) indicate 
that in this specific case, X\ is allowed for both spin- 
conserving and spin-flip transitions, while X4 is only al- 
lowed for spin-conserving transition, originate from the z 
component of the spin-orbit interaction [oc (VV x p) z cr z ] 
that does not flip spin. 

Next we extend the analysis to arbitrary spin orien- 
tation, which leads to the anisotropy of spin relaxation 
processes and enables a direct comparison to a wide range 
of spin injection experiments. 

The spin orientation (n) is defined in terms of polar 
(9) and azimuthal angles (<fi) with respect to the +z and 
+x directions. The new spin states relate to the original 
ones by an 'active' rotation matrix in spin sub-space, 



exp 



cos - 
2 



sin- e 
2 



sin- e 
2 



cos ■ 



(A18) 



where w = n x z/|n x z| is the unit vector along the 
rotation axis. The new spin states follow 



cos-lk^tz) 



sin-e^k,,^), 



(A19) 



sin -e-*\lL L ,1tz) + cos -|k,, ^ 2 ), (A20) 



while the new scattering matrices from k^ to k^ t are 
readily obtained by applying the rotation operator of 
Eq. (|AT8|) on the matrices Eqs. (lAl~6|) and (lAlT)) . The 
new spin-flip matrix elements are 



(kLtAn\Hx?\kL,1tn) 



(1 + i) cos 2 - + (1 - i) sin 2 -e 



= ishx6e i *D x ,„. 



k Lt An\H xb \k L A n ) (A21) 

D Xl , s , 

k Lt A n \H xb \k L ,U) (A22) 



Summing the square amplitudes of the two branches 
leads to Eq. (jHJ in the paper. Table IIIII lists the rela- 
tive magnitudes of the squared spin-flip matrix elements 
for n along several inequivalent high-symmetry directions 
of the crystal. 

The matrix elements are determined by the rele- 
vant directions of the spin orientation and the valley- 
to-valley configurations. For configurations other than 



11 



Lin O ini, the matrix elements could be obtained 
from Eq. by coordinate transformations. If we rewrite 
Eq. ([8|) in the form of the projections of n on x, y, z axes 
as 



\(k Lt A\H X] \k L , 

i=l,2 

2^ iS (l + z 2 + 2^), if j 

,2 

X 4 ,s 



(A23) 



2D 2 y ; „(l - z 2 ), if J = 4, 



then in other valley-to- valley configurations, for example, 
Lyii -H- L^j, the matrix elements are just interchange 
{£, z} of Eq. (IA23[) into {i, — z, y}. Results of all pos- 
sible configurations are listed in Table Q] of the paper. 



TABLE V: Parameters of the L point Hamiltonian [Eqs. (JTH}- 
(I22p ] for bulk germanium following a spin-dependent empiri- 
cal pseudopotential model, mo denotes the free electron mass. 





2.2 


eV 


P 


9 


eV-A 


ml 


0.17m 


Eg iV 


3.3 


eV 


Pi 


7 


eVA 


m*i 


1.60mo 


A L 


0.027 


eV 


P 2 


1.8 


eVA 


mt i3 


1.2mo 


A 3 


0.022 


eV 


a 


0.04 


eV-A 


m l,3 


1.7m 


A 3 , 


0.1 


eV 








m *t,3' 

m *i,y 


-0.16mo 
1.9m 



order terms of intravalley scattering between |ki = K 
q/2, "ft") and |k 2 = K — q/2, JJ.) are found to be 



Appendix B: Derivation of the spin-dependent 
i-point Hamiltonian and the eigenstates 

We use the method of invariants to derive the Hamilto- 
nian. The general procedure is: 1. Figuring out the two 
IRs of the coupling matrix; 2. Decomposing the direct 
product of these two IRs into a sum of IR(s); 3. Accord- 
ing to this decomposition, associating invariant compo- 
nents and matrices to construct the Hamiltonian. These 
invariant components and matrices are obtained by ap- 
plying the symmetry operators on the components of the 
perturbation and the chosen basis states, respectively. 
Associating the invariants to IRs is then carried by exam- 
ining the resulting transformation. Table IIVI lists these 
invariant components and matrices of the Lm point. 
From this table, the constructions of in Eqs. (fT5)l - 
(f2"2")l are straightforward. Parameter constants of these 
Hamiltonian blocks are provided in Table [V] 



TABLE IV: Relevant invariant components and matrices of 
the L m-point. 



IRs 


Invariant components (111) 


Invariant matrices 


L x 


& i ~\(k u + k v ) + k w 


1, / 


L 2 


O w 


Pv 


L 3 


{<T U , &v} 


{[0, l],[l,Q]}, {p x ,p z } 


L' 3 


{-k v ,k u }, {(k x a) u , (k x cr) v } 


{[0, i],[i,o]} 



(£t)®2 4 +,A + _ A +,A £ 



02 



2£U+' A £ 



CIA, 



k fr 



K, it 
(CI) 



Lower order terms (in q) vanish due to the celebrated 
Elliott- Yafet cancelation,^^ In this notation, q® 2 -£® 2 
denotes the scalar product of two second-rank tensors 
(each formed by a dyadic product of the vector with it- 
self). C is the derivative in k-space with its components 
defined by 



£;|k, s) = lim 

<5k^0 



Ik + Sk^s) 



Sh 



(C2) 



In connection with the L point Hamiltonian, C operates 
on the eigenvectors [C 7 (k, s)] and the envelope phase of 
the wavefunction [exp (ik ■ r)]. The electron-phonon in- 
teraction in Eq. (jClj) is given by 



q 



W4 



(C3) 



where the + and — signs denote, respectively, the in- 
phase and out-of-phase motion of atoms in the unit-cell. 
For scattering with long-wavelength acoustic phonon 
modes (A is TA or LA), the out-of-phase polarization vec- 
tor (£q ,A ) is linear in q while the in-phase vector (£+ ,A ) 
has a zeroth-order dependence (e.g., qi/q terms). It is the 
opposite case for scattering with long-wavelength optical 
phonon modes (A is TO or LO). These wavevector depen- 
dencies are taken into account in finding the power-law 
order of the intravalley spin-flip matrix element. Denot- 
ing the atoms positions in the unit cell by ta and tb 
with respect to the cell's origin, the potential in Eq. (|C3[) 
reads 



Appendix C: Wavevector Order Analysis of 
Intravalley Spin-Flip Transitions 

The theory for intravalley spin flips in Ge and Si share 
similar features. In Ref. |13J, we have analyzed the case 
of Si. Here we summarize the important findings and 
discuss the difference for the case of Ge. By invoking 
space inversion and time reversal symmetries the leading 



V±(r)=V at (r - T A ) ± V at (r - r B ), 

V at (r) = Vat (r)I + [ W a t (r) x p] • a , (C4) 

where the spin-orbit interaction is included in the atomic 
potential. 

In the next step of the analysis, we expand the states 
around the valley center. The bra and ket states in 
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Eq. (jCll) are taken at the average of ki = K + q/2 and 
k 2 = K — q/2. We expand this averaged state around 
the valley center position (Ko), 

|K, s) = |K , s) + K • £|K , s) + G(X 2 ) , (C5) 

where K is measured with respect to Ko. Substitut- 
ing this expansion in Eq. (|Clj) . one can identify which 
terms vanish. This identification is carried straightfor- 
wardly using the transformation properties of C, W± 
and |Koj s) under space inversion and time reversal sym- 
metries. 

The difference between the analysis of Si and Ge stems 
from the position of the valley center. The valley cen- 
ter in Ge is at the zone-edge L point, and in Si it is 



inside the Brillouin zone (0.15 x 2ir/a away from the 
X point along the A-axis). Since Ko and — Ko are the 
same point in Ge, space inversion operation keeps |K , s) 
invariant in Ge but not in Si. Together with the transfor- 
mations of C, W±, and |Ko, s) one can readily identify 
the dominant contributions to intravalley spin-flip ma- 
trix elements. Scattering with long-wavelength acoustic 
phonons is led by qtq m products in Si and by Kiq m q n 
products in Ge. For scattering with long-wavelength op- 
tical phonons, the leading terms in Si are linear in q, 
and in Ge they include Kiq m products. Finding the ex- 
act products, their coefficients and deformation potential 
constants requires a combination of kp, rigid-ion and 
group theories.— 
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